Level Splitting in Association with the Multiphoton 
Bloch-Siegert Shift 



P L Hagelstein 1 , I U Chaudhary 2 

1 Research Laboratory of Electronics, Massachusetts Institute of Technology, 
Cambridge, MA 02139,USA 

E-mail: plh@mit . edu 

2 Research Laboratory of Electronics, Massachusetts Institute of Technology, 
Cambridge, MA 02139,USA 

E-mail: irf anc@mit . edu 

Abstract. We present a unitary equivalent spin-boson Hamiltonian in which terms 
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order degenerate perturbation theory is used to develop approximate results in the 
case of moderate coupling for the level splitting. 



PACS numbers: 32.60.+i,32.80.Bx,32.80.Rm,32.80.Wr 



Submitted to: J. Phys. B: At. Mol. Opt. Phys. 



Level Splitting in Association with the Multiphoton Bloch-Siegert Shift 2 
1. Introduction 

The dynamics of a two-level with sinusoidal coupling has been of interest since the 
time of Bloch and Siegert (TJ [2]. The (closely related) basic model for a two-level 
system coupled to a simple harmonic oscillator was considered by Cohen- Tannoudji 
et al [3j. The coupling in these models produces an increase in the two- level system 
transition energy (sometimes termed the Bloch-Siegert shift). As the coupling strength 
is increased, the levels shift relative to one another, producing both level crossings and 
level anticrossings. Level crossings occur when the dressed two-level transition energy 
matches an even number of oscillator quanta (in which case the parity of the states 
are mismatched, so no mixing occurs). Level anticrossings occur when the dressed two- 
level transition energy is resonant with an odd number of oscillator quanta, with the 
magnitude of the splitting indicative of the ability of the coupled system to convert 
energy between the two different degrees of freedom. 

These models were studied initially in the context of spin dynamics in a magnetic 
field [TJ [2] , but they also appear in other applications. The coupling between atoms and 
an electromagnetic field can in some cases be described by these models, in which case 
the resonances mentioned above correspond to multiphoton interactions. Multiphoton 
resonances in which a substantial number of photons are exchanged have become 
experimentally accessible recently pi]. In part because of this there has been renewed 
interest in the multiphoton regime [HIE]. 

We have found a unitary transformation which produces a rotated version of the 
problem which appears to provide a clean separation between terms which produce 
most of the Bloch-Siegert shift, and terms which produce the level splitting at the 
anticrossings. This is interesting because it allows us to develop estimates when the 
coupling is moderately strong for both the shift and the splittings using conventional 
methods on the rotated problem. In essence, we are able to capture most of the level 
splitting in the multiphoton regime in terms of first-order coupling in the context of 
degenerate perturbation theory. This provides a new way to look at the problem which 
may be useful. 
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The Hamiltonian for the coupled two-level system and oscillator of interest (the spin- 
boson Hamiltonian) can be written as 
AE 



H 



-u. 



+ hu a)a + U(a^ + a)a x 



(1) 



where the &i are the Pauli matrices. Since we are interested in the multiphoton regime, 
we assume that the background excitation of the photon field is large: 

AE > hcoo , n > 1 

Rotations are often used to simplify Hamiltonians [TJ; however, in this case our 
rotation will make the problem more complicated mathematically but perhaps simpler 
functionally as outlined above. We consider the unitary equivalent Hamiltonian 

H' = U ] HU (2) 

where 
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The rotated Hamiltonian, H' can be broken up into an unperturbed part (Ho), and 
pieces which will be considered as perturbations (V and W): 

H' = H + V + W (3) 



where 



Hn 



V 



AE 2 + 4U 2 (a + tf) 2 Y + hu tfa 



(4) 



u 

~AE 



1 + 



2U(a + a 
AE 



(a - a f ) 



+ (a — a*) 



U 
~AE 



1 + 



2U{a + qt) 
AE 



(5) 



W = hujo < 
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In the multiphoton regime of interest here (n ^> 1 and AE ^> ^o); the last term 
produced by the rotation, W, is small; we will therefore neglect it in what follows. 
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3. Eigenvalues of H 

Consider first the eigenvalue equation of the unperturbed Hamiltonian H in the rotated 
frame 

Eil> = (7) 
Separation of variables allows us to develop solutions of the form 

■0 = u <S> \s, m) (8) 



where u satisfies 



huo\ htoo 



dy 2 



u(y) + mJAE 2 + 8V 2 y 2 u(y) (9) 



Both from numerical calculations and the WKB approximation we have found that the 
energy eigenvalues are given approximately by 

E n , m (g) = AE(g)m + hu n (10) 

The WKB approximation can be used to develop a useful analytic approximation to the 
dressed two- level transition energy AE(g), which we may write as 



7T J-y/i V e - y 2 
with e = 2n + 1. The dimensionless coupling constant g is 

' ^ U i§ ^ 

In the limit of large n and Huq <^ AE where this is valid, the rotated system governed 
by Ho alone behaves like a dressed two-level system (with increased transition energy) 
and an unperturbed oscillator. 

The condition for Bloch-Siegert resonances can be written as 

AE(g) = (2k+l)hu) (13) 

Level crossings occur in the modified version of the problem described by the 
unperturbed Hamiltonian Hq in the rotated frame. As level anticrossing occur in 
the original spin-boson model at these resonances, the coupling that is responsible for 
the level anticrossing has been eliminated in H . This is an interesting and perhaps 
unexpected feature of this rotation. 
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Near a resonance, we can use a two-level description to account for the level splittings. 

\ I E (g) v W co \ 
v E 1 (g) l\c x ) 

Two levels with energies Eq and E\ that depend on the dimensionless coupling strength 
g cross, and couple to each other with an interaction v which we assume to be constant 
in the vicinity of the resonance. At resonance {go), the two levels in this simplified 
model are degenerate 

E (g ) = Ex (go) (15) 

The splitting between the two levels at this point is twice the magnitude of the 
interaction 

AE min = E+(g )-E4go) = 2\v\ (16) 

The level splittings in the case of weak coupling have been known for some time 
[2], as mentioned above. Shirley's results written in our notation are 

2fc+l //\E\ 2k 

= K#T77^ h- AS (17) 



2 2k - 1 (k\) 2 \huo) 

We have plotted results from the direct numerical solution of the spin-boson Hamiltonian 
[Equation (pQ)], and also for the this weak coupling result in figure [H When the 
dimensionless coupling constant g is small the results match well; when the coupling 
gets stronger, we see (as expected) that perturbation theory begins to break down. 

5. Level splitting in the rotated Hamiltonian 

The dressed transition energy of the two-level system is described reasonably well 
through the unperturbed part Ho of the rotated Hamiltonian, but no level splittings 
occur in the eigenvalues of Hq. Hence, all of the splitting must be due to the terms we 
have considered to be perturbations. In this section, our goal is to apply degenerate 
perturbation theory in the rotated frame to see whether the larger of the perturbation 
terms V can account for the level splitting. 

To calculate the level splitting in the vicinity of an anticrossing, we need to compute 
the eigenkets ^Vm of H Q where 

Ho1pn,m E n ^ m ll) nm (18) 
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Figure 1. Energy level splitting on resonance as a function of the dimensionless 
coupling strength for AE = HJvuiq at large n. Exact numerical results: full circles; 
literature results: open circles. The odd integers label the Bloch-Siegert resonance 
2k+ 1. 



This can be done numerically, or by using the WKB approximation (which we have found 
to be effective for such problems). Near the (2k + l)th resonance, the level anticrossing 
is well- described by a simple two-level approximation 

E t g \ ( C ° ] _ [ E H)m (g) (lpn,m\V\lp n+ 2k+l,m~l) ] | C ° ] 

\ C\ J \ (lpn+2k+l, m -l\V\'>Pn,m} E n+2 k+l,rri~l(g) j \ ° l j 

The energy splitting at resonance is 

AE min = E + (g ) - E^(g ) = 2|(^ nim |V r |^ n+2 fc+i, m -i)| (20) 

In figure [2] we have plotted level splittings taken from a direct numerical solution of the 
original spin-boson Hamiltonian [Equation (pQ)] and also from first-order degenerated 
perturbation theory as discussed here (we used numerical solutions for the eigenfunctions 
ipn,m for this result). We can see from figure [2] that the exact numerical results for the 
level splitting of the unrotated Hamiltonian match very well the results obtained by 
using degenerate perturbation theory on the rotated Hamiltonian. Minor deviations 
occur at the larger g values which we attribute to the omission of higher-order terms in 
the degenerate perturbation theory. 
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Figure 2. Energy level splitting on resonance as a function of the dimensionless 
coupling strength for AE = llhuio at large n. Exact numerical results: full circles; 
first-order splitting from (|20p : open circles. 

6. Conclusion 

We have found a useful unitary transformation that produces a rotated Hamiltonian for 
the spin-boson problem in the multiphoton regime that has interesting properties. The 
rotated Hamiltonian is more complicated mathematically than the initial spin-boson 
Hamiltonian, but appears to be simpler in terms of functionality. One part of the 
rotated Hamiltonian is identified as an unperturbed Hamiltonian (Ho) which appears 
to describe the coupled systems reasonably well away from the level anticrossings. 
This part of the problem is useful for developing estimates of the Bloch-Siegert shift. 
Another part of the rotated Hamiltonian (V) is identified as a perturbation which is 
responsible for most of the coupling which occurs at the anti-crossing. Used with first- 
order degenerate perturbation theory, this term provides a reasonable approximation for 
the level splittings at the Bloch-Siegert resonances. Finally, there is present an additional 
term (W) in the rotated Hamiltonian which is small (so that we have neglected it in 
our discussion here), but which can provide a minor correction to the dressed two- level 
system energies. 
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